A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



A complete sample of GHz-Peaked-Spectrum radio sources 
and its radio properties 



00 



oo 



> 

(N 
(N 

m 
o 

oo 

ON 

Oh- 
I ■ 

o 



C. Stanghellinii ^ CP. O'Dea^, D. Dallacasa'^, S.A. Baum^, R. Fanti^ and C. Fanti 

Istituto di Radioastronomia del CNR, CP. 141, 1-96017 Noto SR, Italy 
Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD, 21218 
Dipartimento di Fisica, Universita degli Studi, via Irnerio 46, 1-40126 Bologna, Italy. 
Istituto di Radioastronomia del CNR, Via Gobetti 101, 1-40100 Bologna, Italy 
^ visitor at the Space Telescope Science Institute, Baltimore 

Submitted 18/12/97 



3,4 



Abstract. We define a complete sample of thirty-three 
GHz-Peaked-Spectrum (GPS) radio sources based on their 
spectral properties. We present measurements of the radio 
spectra and polarization of the complete sample and a list 
of additional GPS sources which fail one or more criteria 
to be included in the complete sample. 

The majority of the data have been obtained from 
quasi-simultaneous multi-frequency observations at the 
Very Large Array (VLA) during 3 observing sessions. Low 
frequency data from the Westerbork Synthesis Radio Tele- 
scope (WSRT) and from the literature have been com- 
bined with the VLA data in order to better define the 
spectral shape. 

The objects presented here show a rather wide range 
of spectral indices at high and low frequencies, including 
a few cases where the spectral index below the turnover is 
close to the theoretical value of 2.5 typical of self-absorbed 
incoherent synchrotron emission. Faint and diffuse ex- 
tended emission is found in about 10% of the sources. 

In the majority of the GPS sources, the fractional po- 
larization is found to be very low, consistept with the 
residual instrumental polarization of 0.3 %. EJ 



Key words: galaxies: active ■ 
continuum: galaxies 



quasars: general — radio 



1. Introduction 

The GHz-peaked-spectrum (GPS) radio sources are char- 
acterized by a simple convex spectrum which peaks in 
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a range of about a decade around 1 GHz. General dis- 
cussions on the general properties of these objects are 
given by O'Dea et al. (1991) O'Dea and Baum (1997), 
and O'Dea (1998) where an exhaustive bibliography can 
also be found. 

Common characteristics of the bright sample of GPS 
radio sources are: small size (^ 1 kpc), high radio lu- 
minosity, low fractional polarization, and apparently low 
variability. They are a mixed class of quasars and galax- 
ies. Galaxies tend to be or brighter and at redshifts 
0.1 ^ ^ ^ 1 (O'Dea et al. 1996) while quasars are often 
found at very large redshift 1 ^ z ^ 4, (O'Dea 1990). 

Currently, there are two main competing hypotheses to 
explain the origin of GPS radio sources and their possible 
evolution. 

In the "young source" scenario, first suggested by 
Philhps and Mutel (1982), GPS radio sources with com- 
pact double (CD) morphology or compact symmetric mor- 
phology (Compact Symmetric Objects: CSO) are classical 
double radio sources at the very first stage of their lives. 

In the "frustration" scenario, GPS radio sources will 
never become as large as the classical doubles since they 
are confined to the sub-kpc scale by a dense and turbulent 
ambient medium (e.g., O'Dea et al. 1991, Carvalho 1994, 
1998). However, it is unclear whether the gas density and 
its distribution in the nuclear region is sufficient to confine 
the radio source for times of the order of 10*" years (see 
e.g., O'Dea 1998). 

Recent results support the young source hypothesis. 
Fanti et al. (1995) presented a model for the evolution of 
the galactic-scale Compact Steep Spectrum (CSS) sources 
into large scale classical doubles. They argued that the 
objects showing symmetric morphology are probably not 
confined by a dense and clumpy medium (see also De 
Young 1993). They also suggest that the typical CSS 
source age is of the order of 10^ years and that the ra- 
dio source luminosity decreases by an order of magnitude 
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as the size of the radio source grows from a few kpc to 
hundreds of kpc (sec also Begelman 1996). The same con- 
clusions are reached by Readhead et al. (1996a,b) in their 
study of a few CSO's. O'Dea and Baum (1997), combining 
the complete sample of GPS presented here with the CSS 
sources by Fanti et al. (1990) reach similar conclusions. 
They further note note that the luminosity evolution re- 
quired makes it likely that GPS and GSS radio sources will 
evolve into objects less powerful than the most powerful 
classical doubles. 

The detection of arc-second scale faint extended emis- 
sion around ~ 10% of GPS sources (Baum et al. 1990, 
Stanghellini et al. 1990) motivated Baum et al. (1990) to 
suggest that nuclear activity is recurrent in these sources. 
In this hypothesis, we see the relic of a previous epoch of 
activity as faint diffuse emission surrounding the current 
young nuclear source. 

In order to achieve a deeper understanding of the GPS 
radio sources we pursued the following project. 

(1) We created a complete sample of GPS radio sources 
by means of a preliminary bibliographic research, checked 
with subsequent new multi-frequency data from the Very 
Large Array (VLA) and the Westerbork Synthesis Radio 
Telescope (WSRT). 

(2) We determined the properties of the radio spec- 
trum and the polarization. 

(3) We obtained optical imaging to determine the host 
galaxy properties. 

(4) We obtained VLBI observations to study the mil- 
liarcsecond radio morphology. 

In this paper we present the selected sample and the 
results from the VLA and WSRT radio observations. The 
optical properties of GPS radio sources are discussed 
by O'Dea, Baum and Morris (1990), Stanghellini et al. 
(1993), and O'Dea et al. (1996). The milliarcsecond mor- 
phology is presented in Stanghellini et al. (1997). Con- 
straints on radio source evolution based on our observa- 
tions of the complete GPS sample are discussed by O'Dea 
& Baum (1997). 

Ho=100 km sec~^ Mpc~^, and qo=0.5 have been used 
in this paper. 

2. The Complete Sample 

Research on GPS radio sources has increased in the past 
few years and a heterogeneous list of objects has been 
accumulated. The list presented by O'Dea et al. (1991) 
consists of ~ 100 radio sources, while Dallacasa and 
Stanghellini (1990) collected a larger list of GPS and CSS 
radio sources and candidates. To enable statistical analysis 
of the properties of these objects we selected a complete 
sample of bright GPS radio sources. 
The selection criteria were: 

- declination 6 > —25° 

- galactic latitude |Z| > 10° 

- flux density at 5 GHz Sqghz > ^Jy 



- turnover frequency between 0.4 and 6 GHz 

- spectral index athin > 0.5 (5*,^ oc i/~") in the high 
frequency, likely optically thin, part of the spectrum. 

We started by selecting GPS candidates from the 1 Jy 
catalog of Kiihr et al. (1981). We cleaned this first "dirty" 
sample (Stanghellini et al. 1990) using our multi-frequency 
observations from the VLA and the WSRT presented here, 
supplemented with data from the literature. 

The final complete sample consists of the 33 objects 
listed in Table 1. This is the first complete sample of bright 
GPS objects. The new radio data permit the proper classi- 
fication and the improvement in the estimate of important 
parameters such as spectral indices, turnover frequencies, 
and polarization properties. 

We note the following regarding the sample and pos- 
sible selection effects. 

Since variability may influence the spectral shape, the 
fact that GPS sources are thought to not be significantly 
variable may be just due to selection effects. Even if the 
spectral shape is constant, if the flux density varies sig- 
nificantly, the source would probably not be recognized as 
GPS using data from the literature spanning several years. 
Therefore it is important to have simultaneous multi- 
frequency observations to build a sample unbiased with 
respect to variability (see also sect. 4.4). 

The turnover frequency and the spectral indices at low 
and high frequencies are difficult to derive when few data 
are available and/or the spectrum bends continuously. For 
this reason we included objects with a turnover frequency 
in a range a little larger than that of the canonical decade 
around 1 GHz. 

The selection is based on the observed turnover fre- 
quency, hence we select objects with different intrinsic 
turnover frequencies at low and high redshifts. A selection 
based on the intrinsic turnover frequency is not currently 
possible due to the lack of complete redshift information 
on large samples of radio sources. 

We note that the sample is unbiased regarding the op- 
tical identification. 



3. The Observations 

The flux densities obtained in the different observing ses- 
sions are presented in Tables 4 and 5. Selected data from 
the literature have been added (with the relevant refer- 
ences) when considered compatible with our observations 
and useful to define the spectral shape. 

During the VLA observing sessions an additional 20 
GPS candidate sources were observed (Table. 2). These 
sources fail one or more criteria of the complete sample. 
Nevertheless, the majority of them show a clear GPS-like 
spectrum. 
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Table 1. Complete sample: column 1 to 10, source name, redshift, optical identification, magnitude and filter/color, flux density 
at 5 GHz, observed and rest frame peak frequency, references for the optical information 1): O'Dea et al. 1991; 2) de Vries et 
al. 1995; 3) Stanghellini et al. 1993; 4) Snellen et al. 1996; 5) Stanghellini et al. 1997; 6) de Vries, private communication; 7) 
White et al. 1993; 8) Hewitt and Burbidge 1987; 9) Hunter et al. 1993; 10) Heckman et al. 1994; 11) Stickel and Kuhr 1993a; 12) 
O'Dea ot al. 1996; 13) Stickel and Kuhr 1993b; 14) Junkkarinen et al. 1991. Last column is a note indicating the radio source 
with extended emission. 



name 


z 


id. 


m 


SsGifz 












ref. 


note 


0019-000 


0.305 


G 


18.4r 


1.10 


0.8 


1.0 


-1.0 


1.2 


3.47 


3,4 




0108+388 


0.669 


G 


22. Or 


1.29 


3.9 


6.5 


-2.1 


0.9 


1.33 


1,3 


e 


0237-233 


2.223 


Q 


16.6V 


3.34 


1.0 


3.2 


-1.9 


0.7 


7.05 


1 




0248+430 


1.316 


Q 


15.5V 


1.24 


5.2 


12.1 


-0.4 


0.6 


1.27 


1 


e? 


0316+161 




G 


23.6R 


2.91 


0.8 


1.2 


-0.7 


1.1 


9.55 


3 




0428+205 


0.219 


G 


19.3R 


2.30 


1.0 


1.3 


-0.6 


0.6 


4.02 


1,3 




0457+024 


2.384 


Q 


19.4V 


1.57 


1.9 


6.3 


-3.5 


0.6 


1.89 


1 




0500+019 


0.583 


G 


21.0i 


1.89 


2.0 


3.2 


-1.6 


0.9 


2.51 


2 




0710+439 


0.518 


G 


18.5i 


1.67 


1.9 


2.8 


-1.3 


0.6 


2.09 


1,3 




0738+313 


0.631 


Q 


16.1V 


3.66 


5.1 


8.3 


-0.7 


0.8 


3.82 


1 


e 


0742+103 




S? 


23r 


3.46 


2.8 


5.7 


-0.7 


0.7 


4.12 


1,5 




0743-006 


0.994 


Q 


17.5V 


2.05 


6.0 


12.0 


-0.6 


0.8 


2.12 


1,5 




0941-080 


0.228 


G 


17.9r 


1.10 


0.5 


0.6 


-0.3 


1.0 


3.40 


3,6 


e? 


1031+567 


0.459 


G 


20.2r 


1.27 


1.3 


1.8 


-0.7 


0.8 


1.87 


1,3 




1117+146 


0.362 


G 


20.1R 


1.00 


0.5 


0.7 


-0.4 


0.8 


3.89 


1,3 




1127-145 


1.187 


Q 


16.9V 


3.82 


1.0 


2.2 


-0.8 


0.6 


5.80 


8 




1143-245 


1.95 


Q 


18.5V 


1.40 


2.0 


5.9 


-2.5 


0.7 


1.69 


1 




1245-197 


1.273 


Q 


20.5 V 


2.31 


0.5 


1.2 


-0.7 


0.9 


8.69 


1 




1323+321 


0.369 


G 


19.2r 


2.35 


0.5 


0.7 


-0.7 


0.6 


7.03 


1,3 




1345+125 


0.122 


G 


15. 5r 


3.03 


0.6 


0.6 


-0.9 


0.7 


8.86 


1,3 




1358+624 


0.431 


G 


19.8r 


1.80 


0.5 


0.8 


-1.5 


0.7 


6.56 


1,3 




1404+286 


0.077 


Sy 


14.6r 


2.66 


4.9 


5.3 


-1.5 


1.6 


2.76 


1,3 




1442+101 


3.544 


Q 


17.8V 


1.19 


0.9 


4.1 


-1.3 


0.9 


2.61 


1 




1518+047 


1.296 


Q 


22.6r 


1.05 


0.9 


2.1 


-0.6 


1.3 


4.58 


1,10 




1600+335 




G 


23. 2r 


2.67 


2.6 


3.9 


-0.2 


0.9 


3.06 


1,3 




1607+268 


0.473 


G 


20.4r 


1.71 


1.0 


1.5 


-1.5 


1.2 


5.44 


1,3 




2008-068 




G 


21.3R 


1.34 


1.3 


2.0 


-0.9 


0.8 


2.64 


1,3 




2126-158 


3.270 


Q 


17.3V 


1.17 


3.9 


16.7 


-1.3 


0.5 


1.23 


1 




2128+048 


0.99 


G 


23.3r 


2.02 


0.8 


1.6 


-1.0 


0.8 


4.93 


1,11 




2134+004 


1.936 


Q 


16.8V 


8.50 


5.2 


15.4 


-1.2 


0.7 


8.59 


8 


e 


2210+016 




G 


22.0i 


1.05 


0.4 


0.6 


-0.6 


1.0 


4.51 


2 




2342+821 


0.735 


Q 


20.1r 


1.28 


0.5 


0.8 


-0.7 


0.9 


6.29 


1,3 




2352+495 


0.237 


G 


18.4R 


1.47 


0.7 


0.8 


-0.2 


0.5 


2.93 


1 





3.1. VLA observations 

The VLA data were obtained in a number of observing 
sessions in the A configuration spanning several years. 
Eighteen candidate GPS radio sources were observed on 
30 December 1984. The results from this first set of ob- 
servations were already published by O'Dea et al. (1990), 
Stanghellini et al. (1990), and O'Dea et al. (1991). We 
consider here only the observed sources belonging to the 
complete sample. 

A second set of VLA observations was obtained on 9 
February 1990, and a third session on 8 June 1991. The 
journal of these observations in shown in Table 3. In the 



third session, where the majority of the sources belong- 
ing to the complete sample were observed, we split the 2 
frequencies in each band to obtain the maximum possi- 
ble separation with adequate sensitivity. This was used in 
the previous observing sessions only in the L band, here 
we used it also in the C, X and P bands. This technique 
proved to be very useful, because it permitted a better 
definition of the spectral shape and a more accurate selec- 
tion of the sources peaking close to the limit of 0.4 GHz 
that we adopted for our sample. 

Furthermore, the observations separated in frequency 
permit us to solve the nn ambiguity that occurs in the 
determination of the Faraday rotation measure. 
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Table 2. Additional objects: column 1 to 10, source name, redshift, optical identification, magnitude and filter/color, flux 
density at 5 GHz, observed and rest frame peak frequency, references for the optical information, references for the optical 
information as table 1. Last column is a note indicating the radio source with extended emission. 



name 


z 


id. 


m 






^rest 








ref. 


note 




0.6 


Q 


21. Or 


1.31 


1.4 


2.3 


-0.2 


0.4 


1.90 


4 




0901 -1-1 1 S 


3.61 


o 


19. 5R 


0.837 






-0.6 


0.2 




1 7 






0.599 


Q 


21. Ir 


2.91 


\J.O 


0.5 


-0.3 


0.5 


8.15 


1 3 






0.844 


o 


19.22V 


1.13 


0.3 


0.6 


-0.3 


0.1 


1.21 


g 






2.07 


o 


19. 2r 


1.96 


1.9 


5.74 


-0.4 


0.1 


2.14 


3 9 








o 


±O.U V 


fi "^1 


Q 


32.1 


-1 P, 




7 sn 


]^ 








Q 

o 


23. Ir 


n fi4n 


n 7 


1.4 


-1 7 


1 n 


9 ^7 


1 ^ 










1 7 nv 


U.OU 1 


1 (S 


4.1 


n 7 


n 7 


1.42 










E 




0.216 


0.6 


1.3 


-0.4 


1.3 


1.04 


1 2 




flQI 4+1 1 4 




G 


20. Or 


0.120 


0.3 


0.7 


-0.1 


1.6 


2.30 


3 




1543-1-005 


0.550 


G 


20.0R 


1.23 






0.1 


0.5 




1,10 




1732-1-094 




G 


20.7R 


0.742 


2.3 


4.5 


-1.1 


1.1 


1.35 


1 




2015+657 


2.845 


Q 


19.7R 


0.667 






-0.3 


0.2 




1,11 




2021+614 


0.2266 


G 


17.9R 


2.82 


8.4 


10.2 


-0.2 


0.6 


3.29 


1 




2050+364 


0.354 


G 


21. 2r 


3.28 


2.1 


2.9 


-0.4 


0.9 


5.82 


1,12 




2137+209 


1.576 


Q 


19R 


0.590 






0.3 


0.9 




1,10 




2149+056 


0.740 


G 


20.4R 


0.840 


2.5 


4.3 


-1.5 


0.7 


1.09 


1,13 




2223+210 


1.949 


Q 


18.2V 


1.07 












1,14 


e 


2230+114 


1.037 


Q 


17.3V 


4.12 


0.5 


1.0 


-0.7 


0.4 


8.35 


1 




2337+264 




G 


20.0i 


0.996 


2.7 


5.5 


-0.2 


0.6 


1.12 


2 





Table 3. Journal of the second and third VLA observing ses- 
sions. Only a few sources were observed at 22 GHz 





date 


VGHz 




09 


Feb 


1990 


1.380 


12.5 


09 


Feb 


1990 


1.630 


12.5 


09 


Feb 


1990 


4.815 


12.5 


09 


Feb 


1990 


4.865 


12.5 


09 


Feb 


1990 


8.435 


12.5 


09 


Feb 


1990 


8.485 


12.5 


09 


Feb 


1990 


22.435 


50 


09 


Feb 


1990 


22.485 


50 


08 


Jun 


1991 


0.302 


3 


08 


Jun 


1991 


0.333 


3 


08 


Jun 


1991 


1.335 


12.5 


08 


Jun 


1991 


1.665 


12.5 


08 


Jun 


1991 


4.535 


12.5 


08 


Jun 


1991 


4.985 


12.5 


08 


Jun 


1991 


8.085 


12.5 


08 


Jun 


1991 


8.465 


12.5 



The data reduction has been performed in a uniform 

way for all the VLA data sets. The errors in the flux den- 
sities are dominated by the calibration errors which are 
estimated to be around 3% at L, C, X bands, around 5% 
at P and K bands, respectively. 



At VLA resolution almost all the objects observed 
are dominated by a single point-like component. The side 
lobes have been deconvolved from the images using the 
AIPS implementation of the Clark CLEAN algorithm 
(Clark 1980, Cornwell and Braun 1988). The data have 
been self-calibrated in phases (Schwab 1980, Cornwell and 
Fomalont 1988) with an initial point model, and when 
subsequent iterations of self-cal have been considered nec- 
essary, we used an appropriate number of non negative 
components obtained from the images in an interactive 
process until convergence to an acceptable solution was 
achieved. In most cases an amplitude self-calibration has 
been performed on the data to improve the final images 
and to allow a search for possible extended emission. 



3.2. The WSRT data 

Westerbork Synthesis Radio Telescope (WSRT) filler ob- 
servations on most of the objects of the complete sample 
have been obtained during 1990 at 327 MHz and during 
1991 at 608 MHz and 327 MHz. Each source has been 
observed in several snapshots a few tens of minutes long. 

The data reduction for the data obtained in the 1990 
has been done with the package DWARF. The data of the 
1991 observations have been reduced with AIPS. These 
latter data in a format readable by AIPS did not have 
the redundant baselines, therefore the data were not self- 
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Table 4. Flux densities for the complete sample, see Table 5 Table 4. continued 
for notes and references 



source 






err. 


ret. 


0019 - 000 


0.325 


2.37 


.24 


4 




0.365 


2.87 


0.04 


31 




0.408 


2.99 


0.14 


8 




0.606 


3.44 


0.20 


9 




0.750 


3.34 


0.16 


9 




1.380 


2.88 


0.06 


1 




1.640 


2.59 


0.06 


1 




4.860 


1.10 


0.02 


1 




22.460 


0.175 


.018 


1 


0108 + 38 


0.325 


0.065 


.003 


11 




0.608 


0.075 


.002 


10 




1.380 


0.427 


.01 


2 




1.630 


0.572 


.02 


2 




2.3 


0.837 


0.03 


6 




2.695 


1.00 


0.02 


9 




3.9 


1.26 


0.04 


6 




4.816 


1.30 


0.04 


2 




4.866 


1.29 


0.04 


2 




7.7 


0.96 


0.03 


6 




8.434 


0.938 


0.03 


2 




8.484 


0.933 


0.03 


2 




11.2 


0.66 


0.02 


6 




22.460 


0.384 


0.03 


2 


0237 - 233 


0.365 


2.96 


0.06 


7 




0.408 


3.67 


0.11 


8 




0.468 


4.76 


0.12 


9 




0.635 


6.20 


0.25 


9 




0.960 


7.02 


0.20 


9 




1.380 


6.19 


0.12 


1 




1.640 


5.91 


0.12 


1 




4.860 


3.34 


0.07 


1 




8.870 


2.36 


0.08 


9 




22.460 


1.09 


0.10 


1 


0248 + 430 


0.302 


0.452 


0.015 


3 




0.333 


0.481 


0.015 


3 




1.335 


0.804 


0.03 


3 




1.665 


0.922 


0.03 


3 




4.535 


1.244 


0.04 


3 




4.985 


1.241 


0.04 


3 




8.085 


1.13 


0.03 


3 




8.465 


1.09 


0.03 


3 




43.0 


0.385 


0.01 


14 


0316 + 161 


0.318 


6.80 


0.30 


9 




0.365 


7.55 


0.11 


7 




0.408 


8.06 


0.25 


8 




0.468 


8.91 


0.22 


9 




0.635 


9.64 


0.24 


9 




0.750 


9.54 


0.24 


9 




0.96 


9.40 


0.30 


6 




1.380 


7.83 


0.20 


2 




1.630 


7.24 


0.20 


2 




2.3 


5.34 


0.20 


6 




3.9 


3.61 


0.03 


6 




4.815 


2.96 


0.10 


2 




4.865 


2.91 


0.10 


2 




7.7 


1.79 


0.08 


6 




8.435 


1.70 


0.07 


2 




8.485 


1.69 


0.07 


2 




11.2 


1.16 


0.03 


6 




22.460 


0.570 


0.030 


2 



source 


VGHz 


o 

ajy 


err. 


ret. 


0428 -1- 205 


0.325 


2.63 


0.08 


4 




0.408 


2.92 


0.15 


5 




0.635 


3.96 


0.27 


9 




0.968 


4.07 


0.13 


9 




1.380 


3.70 


0.10 


2 




1.630 


3.62 


0.10 


2 




4.815 


2.32 


0.07 


2 




4.865 


2.30 


0.07 


2 




8.435 


1.71 


0.05 


2 




8.485 


1.71 


0.05 


2 




22.460 


0.96 


0.05 


2 


0457 + 024 


0.333 


0.16 


0.02 


3 




0.408 


0.33 


0.01 


5 




0.96 


1.40 


0.06 


6 




1.335 


1.74 


0.06 


3 




1.665 


1.84 


0.06 


3 




2.3 


1.97 


0.06 


6 




4.535 


1.65 


0.05 


3 




4.985 


1.57 


0.05 


3 




8.085 


1.19 


0.04 


3 




8.465 


1.15 


0.04 


3 


0500 + 019 


0.325 


0.471 


.025 


4 




0.365 


0.628 


0.029 


7 




0.606 


1.31 


0.09 


9 




0.96 


1.81 


0.10 


6 




1.380 


2.23 


0.07 


1 




1.640 


2.31 


0.07 


1 




2.3 


2.26 


0.08 


6 




3.9 


2.15 


0.03 


6 




4.860 


1.89 


0.06 


1 




8.870 


1.35 


0.05 


9 




11.2 


1.23 


0.03 


6 




43.0 


0.33 


.02 


1 A 

14 


0710 -1- 439 


0.325 


0.56 


0.03 


4 




0.365 


0.655 


0.019 


7 




0.408 


0.75 


0.04 


15 




0.608 


1.27 


0.06 


4 




1.380 


1.94 


0.06 


2 




1.630 


1.95 


0.06 


2 




4.816 


1.68 


0.05 


2 




4.866 


1.67 


0.05 


2 




8.434 


1.28 


0.04 


2 




8.484 


1.28 


0.04 


2 




22.460 


0.66 


0.03 


2 


0738 + 313 


0.302 


0.64 


0.04 


3 




0.333 


0.73 


0.04 


3 




0.608 


1.15 


0.04 


4 




1.335 


1.93 


0.06 


3 




1.665 


2.30 


0.08 


3 




4.535 


3.67 


0.12 


3 




4.985 


3.66 


0.12 


3 




8.085 


3.45 


0.11 


3 




8.465 


3.41 


0.11 


3 




43. 


0.943 


0.03 


14 



6 
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Table 4. continued 



source 


VGHz 


O 


err. 


ret. 


0742 + 103 


0.302 


1.23 


0.06 


3 




0.333 


1.32 


0.06 


3 




1.335 


3.31 


0.10 


3 




1.665 


3.69 


0.10 


3 




2.695 


3.92 


0.20 


9 




4.535 


3.61 


0.10 


3 




4.985 


3.46 


0.10 


3 




8.085 


2.76 


0.08 


3 




8.465 


2.69 


0.08 


3 




20.0 


1.44 


0.15 


19 




37.0 


1.00 


0.09 


16 


0743 - 006 


0.302 


0.31 


0.02 


3 




0.333 


0.38 


0.02 


3 




1.335 


0.70 


0.02 


3 




1.665 


0.81 


0.02 


3 




2.700 


1.40 


0.08 


9 




4.535 


1.99 


0.06 


3 




4.985 


2.05 


0.04 


3 




8.085 


1.97 


0.06 


3 




8.465 


1.93 


0.06 


3 




20.0 


0.99 


0.10 


19 


0941 - 080 


0.302 


3.23 


0.10 


3 




0.333 


3.33 


0.10 


3 




0.96 


3.16 


0.10 


6 




1.335 


2.71 


0.09 


3 




1.665 


2.43 


0.08 


3 




2.3 


1.90 


0.06 


6 




3.9 


1.41 


0.05 


6 




4.535 


1.20 


0.04 


3 




4.985 


1.10 


0.04 


3 




7.7 


0.738 


0.02 


6 




8.085 


0.691 


0.023 


3 




8.465 


0.657 


0.022 


3 




11.2 


0.550 


0.02 


6 


1031 + 567 


0.151 


0.78 


0.04 


27 




0.325 


1.36 


0.07 


4 




0.408 


1.50 


0.07 


5 




0.608 


1.59 


0.08 


4 




1.380 


1.78 


0.05 


2 




1.630 


1.78 


0.05 


2 




4.816 


1.28 


0.04 


2 




4.866 


1.27 


0.04 


2 




8.434 


0.82 


0.02 


2 




8.484 


0.81 


0.02 


2 


1117 + 146 


0.302 


3.58 


0.20 


3 




0.333 


3.73 


0.20 


3 




0.365 


3.82 


0.08 


7 




0.96 


2.96 


0.11 


6 




1.335 


2.46 


0.08 


3 




1.665 


2.21 


0.07 


3 




2.3 


1.67 


0.06 


6 




3.9 


1.24 


0.04 


6 




4.535 


1.02 


0.03 


3 




4.985 


1.00 


0.03 


3 




7.7 


0.70 


0.02 


6 




8.085 


0.642 


0.023 


3 




8.465 


0.613 


0.022 


3 




11.2 


0.51 


0.02 


6 



Table 4. continued 



source 




Sjy 


err. 


rof. 


1127- 145 


0.302 


4.17 


0.20 


3 




0.333 


4.51 


0.20 


3 




1.335 


5.57 


0.18 


3 




1.665 


5.44 


0.17 


3 




4.535 


4.01 


0.13 


3 




4.985 


3.82 


0.13 


3 




8.085 


3.10 


0.10 


3 




8.465 


3.03 


0.10 


3 




14.940 


2.09 


0.30 


7 




20.0 


1.77 


0.18 


19 


1 1 43 — 245 


0.302 


0.175 


0.01 


3 




0.333 


0.224 


0.01 


3 




1.335 


1.518 


0.08 


3 




1.665 


1.654 


0.08 


3 




4.535 


1.47 


0.07 


3 




4.985 


1.40 


0.07 


3 




8.085 


1.03 


0.05 


3 




8.465 


0.99 


0.05 


3 


1245 - 197 


0.302 


7.80 


0.40 


3 




0.325 


7.86 


0.30 


4 




0.333 


8.22 


0.40 


3 




0.365 


8.89 


0.11 


7 




0.408 


8.61 


0.16 


8 




1.335 


5.21 


0.15 


3 




1.665 


4.63 


0.15 


3 




4.535 


2.48 


0.08 


3 




4.985 


2.31 


0.07 


3 




8.085 


1.51 


0.05 


3 




8.465 


1.45 


0.05 


3 


1323 -1- 321 


0.151 


3.78 


0.10 


26 




0.318 


6.32 


0.18 


18 




0.408 


7.02 


0.35 


5 




0.608 


6.47 


0.30 


4 




1.380 


4.88 


0.15 


2 




1.630 


4.45 


0.15 


2 




4.816 


2.35 


0.08 


2 




4.866 


2.35 


0.08 


2 




8.435 


1.63 


0.05 


2 




8.485 


1.63 


0.05 


2 


1345 + 125 


0.178 


4.60 


0.46 


17 




0.302 


7.86 


0.30 


3 




0.333 


8.39 


0.30 


3 




0.365 


8.31 


0.17 


7 




0.408 


8.78 


0.27 


8 




1.335 


5.33 


0.15 


3 




1.665 


4.84 


0.15 


3 




4.535 


3.18 


0.10 


3 




4.985 


3.03 


0.10 


3 




8.085 


2.28 


0.07 


3 




8.465 


2.21 


0.07 


3 




20.0 


1.13 


0.11 


19 
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Table 4. continued 



source 


VGHz 


o 

ajy 


err. 


ret. 


1358 + 624 


0.151 


1.71 


0.04 


27 




0.302 


4.79 


0.15 


3 




0.333 


5.24 


0.17 


3 




0.365 


5.56 


0.07 


7 




0.608 


5.90 


0.20 


4 




1.335 


4.49 


0.15 


3 




1.665 


3.88 


0.13 


3 




4.535 


1.94 


0.07 


3 




4.985 


1.80 


0.06 


3 




8.085 


1.25 


0.04 


3 




8.465 


1.20 


0.04 


3 


1404 + 286 


0.327 


0.177 


0.005 


12 




0.365 


0.179 


0.032 


7 




0.610 


0.245 


.020 


11 




1.380 


0.79 


0.03 


2 




1.630 


1.02 


0.03 


2 




2.300 


1.65 


0.05 


6 




2.695 


2.00 


0.06 


23 




3.900 


2.52 


0.07 


6 




4.816 


2.64 


0.08 


2 




4.866 


2.66 


0.08 


2 




8.434 


1.95 


0.06 


2 




8.484 


1.94 


0.06 


2 




11.200 


1.48 


0.05 


6 




22.460 


0.490 


0.025 


2 


1442 + 101 


0.302 


1.61 


0.05 


3 




0.333 


1.82 


0.06 


3 




0.608 


2.24 


0.10 


4 




0.96 


2.53 


0.09 


6 




1.335 


2.50 


0.08 


3 




1.665 


2.35 


0.08 


3 




2.3 


1.83 


0.07 


6 




3.9 


1.34 


0.04 


6 




4.535 


1.30 


0.04 


3 




4.985 


1.19 


0.04 


3 




'y 
I.I 


0.70 


0.02 


6 




8.085 


0.76 


0.03 


3 




8.465 


0.72 


0.03 


3 




11.2 


0. 1()2 


0.02 


6 


1518 + 047 


0.178 


2.20 


0.44 


17 




0.365 


3.404 


0.139 


7 




0.608 


4.05 


0.20 


4 




0.96 


4.72 


0.24 


31 




1.380 


3.98 


0.12 


2 




1.640 


3.53 


0.11 


2 




3.9 


1.46 


0.05 


31 




4.815 


1.06 


0.03 


2 




4.865 


1.05 


0.03 


2 




8.435 


0.478 


0.015 


2 




8.485 


0.480 


0.015 


2 




11.1 


0.305 


0.03 


31 


1600 + 335 


0.151 


2.17 


0.10 


26 




0.302 


2.46 


0.12 


3 




0.333 


2.45 


0.12 


3 




0.365 


2.47 


0.03 


7 




1.335 


3.05 


0.10 


3 




1.665 


3.11 


0.10 


3 




4.535 


2.75 


0.09 


3 




4.985 


2.67 


0.09 


3 




8.085 


2.22 


0.07 


3 




8.465 


2.17 


0.07 


3 



Table 4. continued 



source 


VGHz 




err. 


ref. 


iOU / -|- ZOO 


U.oio 


o no 


U.Uo 


1 Q 




0.365 


2.53 


0.07 


7 




0.408 


2.95 


0.15 


5 




0.635 


4.51 


0.28 


9 




0.750 


5.16 


0.18 


9 




1.380 


4.86 


0.15 


2 




1.630 


4.42 


0.14 


2 




4.816 


1.72 


0.05 


2 




4.866 


1.71 


0.05 


2 




8.434 


0.96 


0.03 


2 




8.484 


0.96 


0.03 


2 




22.460 


0.293 


0.015 


2 


2008 - 068 


0.325 


1.05 


0.20 


4 




0.365 


1.35 


0.04 


7 




0.608 


1.95 


0.20 


4 




1.380 


2.62 


0.08 


2 




1.630 


2.48 


0.08 


2 




4.816 


1.34 


0.04 


2 




4.866 


1.34 


0.04 


2 




8.434 


0.84 


0.05 


2 




8.484 


0.78 


0.07 


2 




22.460 


0.37 


0.02 


2 


2126 - 158 


0.333 


0.097 


0.010 


3 




1.335 


0.559 


0.018 


3 




1.665 


0.680 


0.022 


3 




4.535 


1.18 


0.04 


3 




4.985 


1.17 


0.04 


3 




8.085 


1.03 


0.03 


3 




8.465 


1.01 


0.03 


3 




90.000 


0.31 


0.03 


13 


zizo -\- U4o 


U.oUz 


o.Ui 




Q 

o 




O.OOO 


O.OO 


U.lo 


6 




U.OOO 


O.O 1 


U.Uo 


/ 




U.4Uo 


4.ii 


U.iz 


Q 

a 




U.OOO 


O.iU 


U.zo 


y 




0.95 


4.95 


0.25 


6 




1.335 


4.08 


0.13 


3 




1.665 


3.73 


0.12 


3 




2.650 


3.16 


0.02 


9 




3.9 


2.50 


0.12 


6 




4.535 


2.16 


0.07 


3 




4.985 


2.02 


0.07 


3 




8.085 


1.37 


0.04 


3 




8.465 


1.31 


0.04 


3 




14.940 


0.82 


0.12 


7 




31.400 


0.43 


0.11 


9 
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Table 4. continued 



bU Lll Ijt; 




Q , 

Ojy 


Grr. 




2134 + 004 


0.302 


0.60 


0.03 


3 




0.333 


0.66 


0.03 


3 




1 


^ 9'i 


nil 


O 






4.62 


n 1 5 


Q 

o 






8 fil 


U.OU 


O 






8 


U.OU 


Q 

o 




o.uoo 


7 4n 


n 9'^ 


Q 
O 






7 9^1 


n 9n 

u.zu 


Q 
O 




90 n 


4 QQ 


u. ou 


1 Q 




4'^ n 

4:0. U 


9 Q9 
z .y z 


n 1 9 

U. -LZ 


1 4 




87 n 

O 1 .u 


1 74 


n ^9 


90 
zu 


2210 + 016 


0.302 


4.23 


0.20 


3 




0.333 


4.53 


0.20 


3 




0.365 


4.54 


0.06 


7 




0.408 


4.65 


0.21 


8 




0.96 


3.17 


0.16 


6 




1.335 


2.84 


0.10 


3 




1 fifi'i 

1 .uuo 


9 4"! 


n 08 

u.uo 


Q 
O 




9 ^ 


1 77 


u.uo 


u 




^ Q 


1 9fi 













n (14 

U.U4: 


Q 
O 








u.uo 


Q 
O 




7 7 


n ^^^^ 


u.uo 


D 




o.uoo 


yJ.yjO 


n 09 

u.uz 


O 






n fi9 


n 09 

u.uz 


Q 
O 




11.2 


0.45 


0.02 


6 


9^49 4- 891 


n 1 ^1 


R4 


n 1 n 

U. lU 


9^ 
zo 




n ^09 


77 


U.OU 


Q 
O 






^ on 


n ^n 

U.OU 


A 




u.ooo 


^ Q1 


U.OU 


O 




n fifiq 


^ 70 


U.OU 


4 




1 ^^^i 
1 .ooo 


o.yu 


n 1 

U. lO 


Q 
O 






98 


n 1 1 

U. -L -L 


Q 
O 






1.41 


u.uo 


o 




4 QS^ 


1 9S 


n n4 

u.u^ 


Q 
O 


2352 + 495 


0.151 


1.45 


0.15 


30 




0.325 


2.50 


0.12 


4 




0.408 


2.61 


0.22 


14 




0.608 


2.82 


0.06 


10 




0.966 


2.66 


0.27 


9 




1.380 


2.59 


0.08 


2 




1.630 


2.43 


0.07 


2 




4.816 


1.48 


0.05 


2 




4.866 


1.47 


0.05 


2 




8.434 


1.07 


0.03 


2 




8.484 


1.06 


0.03 


2 




22.460 


0.683 


0.035 


2 



Table 5. Flux densities for the additional objects. References 
to the flux density measures for tables 4 and 5: 1)VLA first 
session 2)VLA second session; 3)VLA third session; 4)WSRT; 
5)Northern Cross; 6)Ratan 600; 7)Douglas et al. 1996; 8)Large 
et al. 1981; 9) Kiihr et al. (1981); 10) Baurn et al 1990; 11) de 
Bruyn private comunication; 12) de Bruyn, 1990; 13) Steppe 
et al. 1988; 14) Chandler 1995; 15) Ficarra et al. 1985; 16) 
Esko private communication ; 17) Gower et al. 1967, Pilking- 
ton et al. 1965; 18) Dennison et al. 1981; 19) Edelson 1987; 
20) Tcrasranta et al 1992; 23) Fiedler et al. 1987; 24) Steppe 
et al. 1992; 25) Baldwin et al. 1985; 26) Hales et al. 1988; 27) 
Hales et al. 1990; 28) Hales et al. 1991; 29) Hales et al. 1993a; 
30) Hales et al. 1993b; 31) Khabrakhmanov private communi- 
cation, notes: the errors are taken as 3% if they were found in 
literature to be less then 3%. 0108+388: the flux density at 608 
MHz and higher frequencies refers to the core only, while the 
flux density measure at 325 MHz likely includes the extended 
structure found close to the main compact component (Baum 
et al. 1990). 0552+398, the flux density at 37 GHz is the av- 
erage of the data taken in 1990 (Terasranta et al. 1992), and 
the error is the r.m.s of the distribution. The flux density at 
90GHz is the average of the data taken in 1990 (Steppe et al. 
1992), and the error is the r.m.s of the distribution. 2134+004 
The flux density at 87 GHz is the average of the data taken 
between 1985 and 1988 (Terasranta et al. 1992), the error is 
the r.m.s of the distribution. 



source 




S,/y 


err. 


ref. 


0026 + 346 


0.365 


1.36 


0.05 


7 




0.408 


1.88 


0.18 


9 




1.380 


1.91 


0.06 


2 




1.630 


1.87 


0.06 


2 




4.815 


1.32 


0.04 


2 




4.865 


1.31 


0.04 


2 




8.435 


1.05 


0.03 


2 




8.485 


1.04 


0.03 


2 




22.460 


0.745 


0.04 


2 


0201 + 113 


0.323 


0.350 


0.01 


21 




1.380 


0.852 


0.03 


2 




1.630 


0.891 


0.03 


2 




4.815 


0.844 


0.03 


2 




4.865 


0.837 


0.03 


2 




8.435 


0.803 


0.02 


2 




8.485 


0.801 


0.02 


2 




43.0 


0.565 


0.02 


32 



calibrated, producing a lower dynamic range compared 
with the capabilities of this instrument. 

3.3. Other data 

A few objects were observed during the spring and the 
summer of 1991 with the Northern Cross, a radio tele- 
scope near Bologna operating at 408 MHz on the princi- 
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Table 5. continued 



source 


VGHz 




err. 


ref. 


0404 -1- 768 


0.151 


7.12 


0.30 


28 




0.302 


8.50 


0.40 


3 




0.333 


8.12 


0.40 


3 




0.608 


7.22 


0.40 


4 




i.OOO 


0. / i 


U.i / 


Q 

o 




l.DDO 


o.lo 


n 1 c 
U.lD 


Q 
O 




4.0o{3 


o.Uo 


U.IU 


Q 




A Qftc; 

4.yoo 


z.y 1 


u.iu 


O 

o 




ID.U 


i.oU 


n HQ 
U.Uo 


ZO 


0440 — 003 


0.302 


1.20 


.06 


3 




0.333 


1.23 


.06 


3 




1.325 


1.15 


.03 
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Table 5. continued 
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Table 5. continued 
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pie of the Mill's cross (Braccesi et al. (1969) and Ficarra 
et al. (1985)). The radio telescope dates back to the six- 
ties, however, recent mechanical and electronics upgrades 
allow completely automatic observations with increased 
sensitivity. The flux density values have been set to the 
scale of Baars (Baars et al. 1977) adopting a flux density 
of 37.7 Jy for the calibrator source 3C380 (Riley 1988). 

A few additional objects have been observed with 
the Russian radio telescope RATAN 600 (Esepkina et al. 
1979), in transit mode, during various sessions in 1993 
and 1994, at 11.2, 7.7, 3.9, 2.3 and 0.96 GHz (Mingaliev, 
private communication). 

4. Results cind discussion 

In this paper we present the observational results and the 
primary qualitative conclusions. In a future paper we will 
present a quantitative discussion of the results. 

J^.l. The radio spectra 

The flux densities are presented in Tables 4 and 5 and 
plotted in Figures 1 to 5. 

In order to extract information from the radio spec- 
tra, we fitted them with a hyperbola (which is a curve 
tending asymptotically to a straight line at the extrema). 
We first fit the spectral indices in the thick and thin part 
of the spectrum independently, and then fixed these two 
values (which correspond to the angular coefficients of the 
asymptotes) and we performed a least square fit, solving 
for the other parameters of the hyperbola. 

We calculated the frequency of the spectral peak (Ta- 
bles 1 and 2) using the fitted curves. We show the distri- 
butions of the observed and rest frame turnover frequency 
for the 33 sources of the complete sample in Figures 6c and 
6d. In Fig. 6a and 6b we show the high frequency (above 
the peak) and low frequency (below the peak) spectral 
index distribution, respectively. 

4.2. The spectral indices and the turnover frequency 

The GPS radio sources are a mixed group of galaxies and 
quasars, with some remarkable differences between the two 
classes. The histogram in Fig. 16a shows that the redshift 
distribution is very different for galaxies and quasars. The 
galaxies have a typical redshift of ~ 0.5, and none has a 
redshift higher than 1. For the galaxies without redshift 
information, we note that only 0316+161 has an optical 
magnitude slightly fainter than the galaxy 2128-1-048 at 
redshift 0.99 (Table. 1). Since these galaxies follow the 
Hubble diagram (O'Dea et al. 1996; Snellen et al. 1996) 
it is unlikely they will be found at a redshift much higher 
than the others in the sample. The quasars are instead 
found at any redshift (we included the galaxy 1404+286 
(OQ208), which has a Seyfert 1 nucleus, in the quasar 
class) with the majority at very high z (see also O'Dea 
1990). 
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Fig. 1. Radio spectra of the radio sources of the complete sample 
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Fig. 2. Radio spectra of the radio sources of the complete sample 
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Fig. 3. Radio spectra of the radio sources of the complete sample 



The high frcqiKnic;y spectral index ranges from 0.5 (the 
limit set in the selection criteria) to 1.3 with the galaxies 
having perhaps slightly steeper values (but the 2 objects 
with the steepest spectral indices are quasars). The low 
frequency spectral index ranges from -0.2 to -2.1 without 
any clear difference between galaxies and quasars, but the 
higher values are biased since in several objects the low 
frequency part of the spectrum is under-sampled and the 
spectral index is calculated close to the turnover frequency 
where it is likely to be flatter. Similar results were; found 
by De Vries et al. (1997) though their poorer frequency 
coverage resulted in a somewhat smaller range in spectral 
index. 



The qiiasars tend to peak at higher frequencies than 
the galaxies in both the rest frame and observed frame and 
some quasars have a turnover frequency in their rest frame 
exceeding 10 GHz. This suggests that, on the assimiption 
that the turnover is caused by synchrotron self-absorption, 
the GPS quasars are more compact than the galaxies. This 
effect has been also found by De Vries et al. (1997) in a 
bright heterogeneous sample and by Snellen (1997) in a 
fainter sample. In addition, VLBI images of several sources 
belonging to the complete sample show that quasars are 
more compact than galaxies, and in general exhibit differ- 
ent morphologies (Stanghellini et al. 1997). 

These results suggest that either GPS galaxies and 
GPS quasars are different types of objects, or that beam- 
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Fig. 4. radio spectra of the additional radio sources 
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Fig. 5. radio spectra of the additional radio sources 



ing of compact components plays a role in the quasars (see 
also O'Dea 1998). 

4-3. Extended emission 

We have detected extended emission (both diffuse and 
compact) close to the compact radio source in some cases 
at 21 cm. In the remaining sources, our upper limits on 
extended emission is typically 1 mJy/beam at 21 cm. 

0248-1-430 (Figure 7) has a compact emitting region 
15 arcsec east of the main component and a hint of weak 
emission 5 arcsec to the south. 0528+134 is resolved, show- 
ing an extension in the NW direction (Figure 8). Murphy 
et al. (1993) present an image of 0738-1-313 at 20 cm show- 



ing 2 emitting regions resembling 2 weak hot-spots and 
lobes on the opposite sides of the dominant component. 
In our image (Figure 9) these 2 weak components are al- 
most completely resolved out and only a hint of emission 
has been detected 30 arcsec north and south of the com- 
pact region. 0941-080 shows a slightly resolved secondary 
component 20 arcsec cast of the main one (Figure 10). 
2134-1-004 has very weak and diffuse emission around the 
strong compact component (Figure 11). 2223-1-210 has a 
secondary component 4 arcsec away from the main one in 
the SW direction in our image at 1.35 GHz (Figure 12); 
the main component itself is resolved in a core-jet struc- 
ture oriented NE with a possible counter jet in the image 
at 5 GHz (Figure 13). 2230-^114 at 1.35 GHz (Figure 14) 
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Fig. 6. Histograms for the complete sample: a) high frequency spectral index distribution; b) low frequency spectral index 
distribution; c) observed turnover frequency distribution; d) rest frame turnover frequency distribution 



shows an elongated structure in the NW-SE direction with 
a hint of emission bending to SW, while in the 4.9 GHz 

image (Figure 15) the elongated structure turns out to be 
a core-jet structure with the possible presence of a counter 
jet. 

Stanghellini et al. (1990) report several cases of ex- 
tended emission around GPS radio sources. Of the objects 
presented here showing extended emission, 0528-1-134 and 
2223-1-210 are not true GPS objects (see also section 4.4 
for a discussion of the case of 0528-1-134). In a couple of 
sources (0248-^430, 0941-080 both belonging to the com- 
plete sample) it is difficult to say whether the secondary 
emission is related to the GPS radio source and further 



observations are probably needed. The extended emission 
found around 0738-1-313, 2134-1-004, and 2230-1-114 (the 
first 2 objects belong to the complete sample) is likely to 
be related to the GPS object. In the complete sample, 
0108-1-388 is known to have extended emission, so there 
are 3 to 5 objects out of 33 with known extended emission 
so far. This percentage of 9 to 15% is slightly smaller but 
consistent with that previously claimed by Stanghellini et 
al. (1990). It is clear that the vast majority (^ 90%) of 
the GPS sources appear to be truly isolated and have no 
emission beyond the kpc scale at the current limits. 
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Fig. 16. Histograms for the complete sample: a) redshift distribution; b) fractional polarization at 1.3 GHz; c) fractional 
polarization at 4.9 GHz; d) fractional polarization at 8.5 GHz 



4-4- Variability 

Waltman et al. (1991) presented monitoring observations 
at 2.7 and 8.1 GHz for several GPS sources covering the 
time range from 1979 to 1988. Some sources as 0237-233, 
1245-197, 1345+125 wore foimd to bo very stable in flux 
density. Others were found to be variable: 0552+398 shows 
a variation of ~ 30 % at 8.1 GHz. 2134+004 has a variabil- 
ity of about 15-20% at 8.1 GHz. 2352+495 has a variabihty 
below 10% at 8.1 GHz. Also 0742+103 is slightly variable. 

Wehrle et al. (1992) also report variability for some 
GPS objects in the time range 1985-1991 at 4.8, 8, and 
14.5 GHz, from the University of Michigan Radio Astron- 



omy Observatory monitoring program for several sources, 
some of which are GPS objects. 0552+398 shows an in- 
crease in flux density exceeding 50% at 8.4 and 14.5 GHz, 
1127-145 shows a quasi periodical variability of approxi- 
mately 1 Jy at all the 3 frequencies. The source 2230+114 
also shows a rather remarkable flux density variability at 
all the 3 frequencies with an amplitude of 0.5-1 Jy and is 
a well known low frequency variable source (Bondi et al. 
1996). The variability of 1404+286 has been discussed by 
Stanghellini et al. (1996). 

In conclusion, we find that some GPS sources (mainly 
quasars) show mild to high flux density variability at cm 
and mm wavelengths. However, without uniform moni- 
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toring of the complete sample it will not be possible to 
determine how common this variability is. We also note a 
couple sources where the spectral shape is variable and 
at some; timcis the spectrum was peaked, and at other 
times it was not, 0528+134, the well known gamma-ray 
source (Mukherjee et al. 1996), and 0201+110. Both of 
these sources show a rather flat spectrum in the VLA ob- 
servations from the second or the third session and they 
would be very easily discarded as GPS radio sources. But 
0528+134 has been included in the class of GPS radio 
sources because of its GPS-like spectrum from the liter- 



ature (O'Dea et al. 1991), and 0201+113 really shows a 
convex spectrum in the data in the first VLA session pub- 
lished in O'Dea et al. (1990). This behavior is not surpris- 
ing in highly variable radio sources as we may well expect 
that the presence of new radio components will change 
the spectral shape. Thus, there are sources which show a 
peaked spectrum only part of the time. This aspect of the 
GPS phenomenon deserves more attention as it could be 
related to the remarkably diflferent properties found be- 
tween GPS galaxies and (some?) GPS quasars. 
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4-5. Polarization 

Due to the low level of observed polarization, the errors are 
dominated by the r.m.s noise on the polarization images 
(typically 0.5 mJy) and by contamination from residual 
unpolarized emission (estimated about 0.2-0.3% of the to- 
tal flux density). The error in the polarized flux may then 
be calculated as ap = -^/O.S^ + (0.003 x SmjyYmJy. We 
considered any measurement of the fractional polarization 
below 0.3% to be an upper limit (Tables 6 and 7). The er- 
rors in the position angle have a systematic contribution 
due to the uncertainty in the determination of the position 
angle of the cahbrator (3C286), assumed to be around 2-3 
degrees. This is the dominant contribution in the position 
angle error for most of the objects with detected polarized 
flux density. 

In Fig. 16 we show the histograms of the fractional of 
polarization (mostly upper limits) for the complete sample 
at 1.3, 4.9 and 8.5 GHz. The fractional polarization is in 
general low at all the frequencies. Only a few quasars have 
a fractional polarization above 1 % at 4.9 or 8.5 GHz. 

When polarized emission has been detected we at- 
tempted a linear fit to the polarization angles versus the 
squared observed wavelength, as is expected from the 
Faraday effect on polarized radiation propagating through 
a magnetized and ionized medium. 



The low level or even the lack of detection of polarized 

flux limited us to only a few sources (all quasars) . The fits 
are generally rather good and are given in Table 6 and 7, 
and in Figures 17 and 18. Due to the better frequency cov- 
erage in the present observations, our estimated rotation 
measures supersede those reported by O'Dea et al. (1990), 
though we cannot rule out that some of the difference is 
due to variability. We find Faraday rotation measures in 
the rest frame above 1000 rad/m^ for 5 quasars of the 
complete sample. We also found a very high value (> 10^ 
rad/m^) for 0552+398 which does not belong to the com- 
plete sample but has a GPS shape. 

Sometimes the frequencies which give a good flt in- 
clude those close to the turnover, and in the case of 
0552+398 are all below the turnover. This implies that 
the region emitting the polarized emission is different from 
that responsible for the optically thick emission or that the 
turnover is not caused by synchrotron self absorption. 

4-6. Summary and conclusions 

We have presented a bright flux-density-limited complete 
sample of 33 GPS radio sources selected on the basis 
of their peaked radio spectra. The sample selection was 
based on observations with the VLA, WSRT, and other 
instruments. Additional GPS sources not belonging to the 
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Fig. 7. 0248 + 430 at 1.35 GHz. The restoring beam is 1.41 x 
1.32 arcsec in P.A. —78°. The r.m.s. noise on the image is 0.2 

mjy. The peak flux is 810 mjy/beam. The contour levels for 
all the images are -3, 3, 6, 12, 25, 50, 100, 200, 500, 1000 x the 
r.m.s noise 
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Fig. 8. 0528 + 134 at 1.66 GHz. The restoring beam is 1.14 x 
1.01 arcsec in P.A. +2°. The r.m.s. noise on the image is 0.25 
mJy. The peak flux is 2088 mJy/beam. 



Fig. 9. 0738 + 313 at 1.33 GHz. The restoring beam is 2 x 2 
arcsec. The r.m.s. noise on the image is 0.3 mJy. The peak flux 
is 1943 mJy. 
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Fig. 10. 0941 - 080 at 1.33 GHz. The restoring beam is 2.44 
X 1.54 arcsec in P.A. +29°. The r.m.s. noise on the image is 
0.5 mJy. The peak flux is 2088 mJy/beam. 



complete sample have also been observed. We present 
our results on the polarization and radio spectrum. We 
found remarkable differences in the properties of quasars 
and galaxies, the latter having lower turnover frequencies, 
mostly undetectable polarization and lower redshifts. 

In the few objects where polarization has been de- 
tected at many frequencies, the Faraday rotation measure 
in the rest frame often exceeds 1000 rad/m^. 

In about 10% of the sources we detect weak diffuse 
extended emission. In the remaining ~ 90% any extended 
emission has a peak surface brightness is less than about 
1 mJy/beam at 21 cm. 
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Fig. 11. 2134 + 004 at 1.33 GHz. The restoring beam is 1.91 
X 1.55 arcsec in P.A. —5°. The r.m.s. noise on the image is 0.4 
mjy. The peak flux is 3249 mjy/beam. 



Fig. 13. 2223 + 210 at 5 GHz. The restoring beam is 0.37 x 
0.35 arcsec in P.A. +5°. The r.m.s. noise on the image is 0.3 
mJy. The peak flux is 1034 mJy/beam. 
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Fig. 12. 2223 -h 210 at 1.33 GHz. The restoring beam is 1.37 
X 1.32 arcsec in P.A. —7°. The r.m.s. noise on the image is 0.3 
mJy. The peak flux is 1663 mJy/beam. 



Fig. 14. 2230 -h 114 at 1.33 GHz. The restoring beam is 1.46 
X 1.36 arcsec in P.A. -|-1°. The r.m.s. noise on the image is 0.8 
mJy. The peak flux is 6444 mJy/beam. 
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Fig. 15. 2230 + 114 at 5 GHz. The restoring beam is 0.38 x 
0.36 arcsec in P.A. —11°. The r.m.s. noise on the image is 0.5 
mjy. The peak flux is 4137 mjy/beam. 



In a following paper we will discuss the implications 

of the properties of the complete sample in the framework 
of the scenarios proposed to explain the existence of the 
GPS radio sources. 
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